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ABSTRACT 

Asymmetric nuclear matter is treated in the formalism of Dirac-Brueckner approach 
with Bonn one-boson-exchange nucleon-nucleon interaction. We extract the symmetry 
energy coefficient at the saturation to be about 31 MeV, which is in good agreement 
with empirical value of 30 ±4 MeV. The symmetry energy is found to increase almost lin- 
early with the density which differs considerably from the results of non-relativistic ap- 
proaches. This finding also supports the linear parameterization of Prakash, Ainsworth 
and Lattimer. We find, furthermore, that the higher-order dependence of the nuclear 
equation of state on the asymmetry parameter is unimportant up to densities relevant 
for neutron stars. The resulting equation of state of neutron-rich matter is used to cal- 
culate the maximum mass of neutron star, and we find it to be about 2.1M & . Possible 
mechanisms for the softening of the equation of state are also discussed. 



1 Introduction 

Although the fact that the equation of state of nuclear matter contains a symmetry energy term 
has been known since the early days of nuclear physics, the experimental and theoretical study 
of the symmetry energy and its density dependence is becoming an increasingly interesting topic, 
mainly because of the recent development of radioactive ion beam facilities that allow one to study 
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the structure and reactions of neutron-rich nuclei |l], |2j, |J| , in which the symmetry energy plays an 
important role. The recognization that the symmetry energy, especially its density dependence, 
has a profound effect on the properties of neutron stars [Q, [|, |6|, [7|, || also makes the experimental 
and theoretical determination of this quantity very relevant and useful. 

Empirically, the symmetry energy coefficient S2 (po) in nuclear matter at the saturation density 
po can be extracted from the systematic study of the masses of atomic nuclei, based on, e.g., 
the liquid droplet model |], [l0| or the macroscopic-microscopic model [|ll|, p^j . This, however, 
determines the symmetry energy only for small asymmetry parameter a (a = (N — Z)/A) and 
for densities around pq. From the experience with symmetric nuclear matter we know that the 
determination of the compression modulus K at po does not uniquely constrain the equation of 
state at high densities. Similarly, the determination of S^po) does not guarantee an unambiguous 
determination of symmetry energy at high densities which is needed for the study of neutron star 
properties. 

The situation changes with the recent advances in the development of various radioactive ion 
beam facilities around the world that will produce nuclei with large neutron excess near and beyond 
the drip-line. The study of the structure of these neutron-rich nuclei allows us to determine the 
symmetry energy for large asymmetry parameter and extract possible higher-order dependence on 
a. Furthermore, the collisions of neutron-rich nuclei at relativistic energies, during which nuclear 
matter with densities up to (2-3) po is created, make it possible to study experimentally the density 
dependence of the symmetry energy 14]. 

On the theoretical side, the symmetry energy has been studied over many years based on various 
models and approaches, which can roughly be divided into phenomenological and microscopic, each 



of which can be subdivided into relativistic and non-relativistic. Hartree-Fock [15] and Thomas- 
Fermi |l(| calculations with Skyrme-type effective nucleon-nucleon interactions have been carried 
out to study various nuclear properties, including symmetry energy. The symmetry energy coef- 
ficient Sq from these calculations ranges from about 27 to 38 MeV, and is in agreement with the 



empirical value of 30 ± 4 MeV [17]. Since the the parameters in the Skyrme forces are adjusted to 
fit nuclear matter and finite nuclei properties, this degree of agreement is fully expected. 

Another phenomenological approach that has been used extensively in the study of nuclear 
properties is quantum hadrodynamics (QHD) which is based on the relativistic field theory |18L 
19]. The symmetry energy in this approach has two contributions [EH, 21]; one comes from the 



'kinetic' energy difference between symmetric and asymmetric matter, and the other arises from 
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the exchange of the rho meson which couples to neutron and proton with opposite sign. The 



symmetry energy in this approach ranges from about 35 to 40 MeV [21, ^2], p3f| , somewhat larger 
than the emprical value of 30 ± 4 MeV. 

On a more microscopic level, it has been the ultimate goal of traditional nuclear physics to de- 
scribe in a consistent way the properties of nuclear matter, finite nuclei, and nuclear reactions from 
realistic nucleon-nucleon interactions that are fitted to the nucleon-nucleon scattering and deuteron 
data. There are basically two approaches to this; namely, variational-type and Brueckner-type. In 
variational calculations, it is well-known that the realistic two-nucleon potential alone, such as 
Argonne U14 (AV14) or Urbana v\4 (UV14), does not provide a satisfactory description of nuclear 
matter properties (24], p5| . A phenomenological three-nucleon interaction has to be introduced, 
whose parameters are adjusted so that the variational calculations with the two- and three-body 
interactions give correct nuclear matter saturation density, binding energy and compression mod- 
ulus. The symmetry energy coefficient obtained in the variational calculations is about 30 MeV 



mm- 

The Brueckner-type approach has both non-relativistic and relativistic versions, known as 
Brueckner-Hartree-Fock (BHF) @ |7| and Dirac-Brueckner-Hartree-Fock (DBHF) @ ||, |(| |l, 
|32|j , respectively. As in the variational calculations, the BHF calculation with realistic two-nucleon 
interactions such as Bonn and Paris potentials does not provide a good description of nuclear mat- 



ter properties [26|. On the other hand, the DBHF approach, owning to the additional density 
dependence introduced through the medium modified Dirac spinors, does provide a very good de- 
scription of nuclear matter properties based entirely on realistic nucleon-nucleon interactions that 
are constrained by two-nucleon data (nucleon-nucleon scattering and deuteron properties). 

In this work, we study systematically the properties of asymmetric nuclear matter in the formal- 
ism of the DBHF approach using the Bonn one-boson-exchange (OBE) potential. We will calculate 
in particular the density dependence of the symmetry energy that is very important for neutron 
star properties and heavy-ion collisions. Our results will be compared to those of the variational 
and BHF calculations, and to various phenomenological parameterizations. We will also apply the 
nuclear equation of state obtained in this study to calculate neutron star properties. In Section 2, 
we review briefly the DBHF approach and its extension to asymmetric nuclear matter. The results 
for symmetric energy will be presented in Section 3, while those for neutron stars in Section 4. A 
short summary is given in Section 5. 
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2 Dirac Brueckner approach for asymmetric nuclear matter 



The essential point of the DBHF approach is the use of the Dirac equation for the description of the 
single-particle motion in the nuclear medium. The Dirac spinor, which enters the evaluation of in- 
medium nucleon-nucleon potential, becomes density dependent. This additional density dependence 
is instructive in reproducing correctly the nuclear matter saturation density and binding energy 



The basic quantity in the DBHF calculation is the G matrix which satisfies the in-medium 
Thompson equation, 
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(m{k)\ Q(k,P) 
\E(k)J 2E(q)-2E{k) 



G(q',q\P,z) = V(q',q)+ / W *,Y"'1L^ G (M^*) W 



where E = y/rh 2 + (P/2 + k) 2 , and rh = m + U$, with m being nucleon mass in free space. For 
asymmetric nuclear matter, the angle-averaged Pauli-blocking operator has to be modified and is 
given by 



Q(k,K) 



1 



if Pn>l 



(1 + (3 n )/2 if -1< (3 n < 1 and P P > 1 
(/?„ + (3 p )/2 if fa < 1 and < (0 n + [3 p )/2 
if (/? n + /3 p )/2 < or P n <-1 



(2) 



where 



fin,p 



K 2 /i + k 2 -kl 



Kk 



(3) 



where kp n and kp v are neutron and proton Fermi momenta, respectively, with kp n > fc_F p - 
From the G matrix we can calculate the single-particle potential 



(4) 



which is usually parameterized in terms of the scalar and vector potentials 

fh{k) . 



S(fc) 



E{k) 



U s (k) + U v (k). 



(5) 



In principle, the scalar and vector potentials are both density and momentum dependent. The 
momentum dependence for momenta below the Fermi momentum is rather weak and smooth fl29|, 
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30, 33]. In a recent work (3J|, we found that the saturation properties and nuclear equation of state 
of symmetric matter with explicit momentum dependence are almost the same as those obtained 
in Ref. [3C] under the assumption of momentum independent mean fields. In this work we thus 
follow Ref. [[33] and neglect the momentum dependence of the scalar and vector potentials. 



We will use the Bonn A potential from Ref. [30]. The parameters in this potential are fitted 
to neutron-proton {np) scattering and deuteron properties. It is well known that at low energies, 
there are both charge independence breaking and charge symmetry breaking 133], namely, np, pp, 
and nn potentials are not completely identical. The effects of those differences on the bulk nuclear 
matter properties are, however, found to be extremely small and can well be neglected 

In the case of asymmetric nuclear matter, the potential energy of a single particle is 

E pot = j- / d 3 k-Z n (k) + / d 3 k-E p (k) , (6) 

fo Fn d?k + / Fp d 3 k V° 2 Jo 2 J 

while the kinetic energy is given by 

1 / r k F n 3 mm*(k) + k 2 f k F P 3 mm*(k) + k 2 \ 

km = jS**k + jZ**k U &® + h E*(k) J' (7) 

The energy per nucleon, or nuclear equation of state, is then given by 

E = E pot + Ef, in — m. (8) 

The asymmetric nuclear matter is characterized by the asymmetry parameter 

a= Pn^Pp^ (g) 
Pn + Pp 

In the left panel of Fig. |l] the neutron and proton scalar and vector potentials are shown as a 
function of the asymmetry parameter, for three different densities. It is seen that the variation of 
the proton and neutron potentials is not symmetric with respect to their common value at a = 0. 
For a nucleon at rest in nuclear matter, we can define its single-particle potential as Us + Uy . We 
find that as the asymmetry parameter increases, the neutron potential becomes less attractive and 
that of proton becomes slightly more attractive. These differences in proton and neutron potentials 



may be observed in heavy-ion collisions with radioactive ion beams [13, 14]. 

In the right panel of Fig. [I], we show the nuclear equation of state for a number of asymmetry 
parameters. We compare our results with those of Ref. [^] obtained in the BHF approach. As is 
well-known, the BHF approach saturates nuclear matter at a much too high density. In the DBHF 
calculation, the nuclear matter saturation properties are better reproduced. The binding energy 
and the saturation density become progressively smaller as a increases. 
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3 Symmetry Energy of nuclear matter 



In this section, we discuss our results for symmetry energy and its density dependence. We introduce 
AE as the energy difference between symmetric and asymmetric nuclear matter, 



AE = E(p,a)-E(p,0), 



(10) 



with E(p,a) defined in Eq. (7). The results are shown in the left panel of Fig. |2| as a function of 
a 2 for six different densities. It is seen that at all densities considered here, AE increases almost 
linearly with a 2 , indicating that the a 4 and higher-order terms are not important. To a good extent 
we can express the equation of state of asymmetric nuclear matter as 



E(p, a) = E(p, 0) + S 2 (p)a 2 + S 4 (p)a 4 . 

The usual symmetry energy S 2 is thus defined as 

^d 2 E(p,a) bin 



similarly, 



da 2 



1 d*E(p,a) bin 



a=0 



(11) 



(12) 



(13) 



a=0 



24 da A 

The density dependences of 52 and S± obtained in our calculation are shown in the right panel 
of Fig. [2] by the solid curve. S 2 increases almost linearly with density. Actually a parameterization 
in terms of {p/po) ' 9 fits the theoretical curve reasonably well, as shown in the figure by dotted 
curve. At nuclear matter saturation density, our calculation gives a symmetry energy coefficient 
5*2 (po) °f about 31 MeV, which is in good agreement with the empirical value of about 30 ± 4 
MeV [O] . The BHF and the variational calculations also reproduce the empirical symmetry energy 



coefficient [25, 26]. The coefficient of a 4 term is very small in the density region considered here. 
This means that the approximation of neglecting this term as adopted in Ref. [24, 
reasonable. 



is quite 



In the left panel of Fig. || we compare our results for the density dependence of the symmetry 
energy with those of Ref. pfj] based on the BHF calculations, and of Ref. [25] based on the varia- 
tional calculation. There are significant differences between the results of these three calculations. 



In relativistic approaches |35|, 36], the symmetry energy increases almost linearly with density, and 
is considerably larger than those in non-relativistic and variational calculations [25, |2(|. The differ- 
ence between the DBHF and BHF is mainly due to the relativitic effects. In the simple mean-field 
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approximation to the Walecka-type model, the symmetry energy has a contribution from the 'ki- 



netic energy' difference, which is inversely proportional to E F = ykp + m* 2 . This contribution is 
thus larger in relativistic approaches because of the dropping nucleon mass. This also accounts for 
part of the difference between our results and that of Ref. p!| which is non-relativistic in nature. 
The remaining difference can be explained by the differences in the nucleon-nucleon potentials used 
in the two calculations. In the variational calculations 25], the major contribution to the 'poten- 
tial' part of the sysmetry energy comes chiefly from the second-order tensor interaction, which is 
progressively blocked with increasing density With the strong p-coupling of the Bonn potential, 
the second-order tensor force is relatively weak, compared with that of Ref. [|25|] , so that this is 
not a large effect in our calculation, where the main contribution to the symmetry energy comes 
from /3-meson exchange. The differences in the symmetry energy in these three calculations will 
have a profound impact on the properties of neutron stars. We hope that future experiments with 
radioactive ion beams will help to shed light on this problem. 

Phenomenologically, Prakash, Ainsworth and Lattimer || proposed the following parameteri- 
zation for the density dependence of the symmetry energy, 

S(u) = (2 2 / 3 - 1)\Ep(u 2 I z - F{u)) + S F(u), (14) 
5 

with F(u) = 2u 2 /(l + u), u, or y/u, where u = p/po and E F is Fermi energy at saturation density 
Pq. In the right panel of Fig. ||, the density dependences of three forms of F(u) are compared with 
our results, and it is seen that the F(u) = u case is very close to our results. 



4 Nucleon star properties 

In this section we apply the nuclear matter equation of state obtained in the DBHF calculation to 
the study of neutron star properties. In a real neutron star, there exists a large fraction of protons 
with electrons and muons maintaining the charge neutrality. We call this star a nucleon star instead 
of neutron star ]3^]. We assume that the star is in chemical equilibrium at zero temperature after 
the neutrinos leave the system. The equilibrium condition requires that the chemical potentials 
should satisfy 

P = Pn ~ Pp = Pe = Pfi, (15) 
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where muons start to contribute when p e > m„. The nucleon chemical potential difference p can 
be calculated once we have the coefficients of symmetry energy from DBHF [pSfl , 

A = 4(1 - 2x) \s 2 (p) + 2S 4 (p)(l - 2x) 2 } , (16) 



where x being the proton fraction. From the charge neutrality p e + Pa = Pp and Eq. fll5|), we have 

1 o 1 / 9 9 \ 3/2 

^ = ^2^e +6 '(/ ie ~ m M)^2 (^e • ( 17 ) 

By solving Eqs. ((l^) and ([1?]) for a given density, we can get the proton fraction in nucleon star 



matter. In Fig. |], the proton fraction and the energy difference (Eq. (10)) are plotted in the left 
panel. Muons start to contribute almost at the nuclear saturation density, po = 0.16 fm -3 , which 
leads to a slight bump in the proton fraction. 

Given the symmetry energy and EOS of symmetric nuclear matter, we can get the energy 
density and pressure and solve the Tolman-Oppenheimer-Volkov (TOV) equations for a spherically 
symmetric nucleon star, 

dM A 2 

= 47rr e 

dr 

dP _ G[e + P][l+ 47rr 3 P] 

dr r{r-2GM) ' [ ' 

In the above equations, the pressure P{e) determines the stiffness of the TOV equation, On the 
right panel of Fig. the nucleon star properties are summarized, where the dash-dotted curves 
give the results without inclusion of muons. The maximum mass nucleon star corresponds to 
M max = 2.1 Mq, R = 10.8 km, and the central density p c = 1.08 fm -3 . This mass, we believe, 
is well determined in our conventional scenario, beginning from the microscopic Bonn two-body 
interactions, and using the many-body theory as described above. The EOS can be substantially 



softened, however, by introduction of kaon condensation [39, 40 1 and this brings the maximum mass 



down to M max ~ 1.5 Mq [41|. Introduction of K~-mesons implies an ~ 50% nucleus components 
of protons at the higher densities, so that "nucleon" star, rather than neutron star, then becomes 
the most appropriate name. A recent compilation by Steve Thorsett quoted by Brown |[7|] shows 
that well-measured neutron star masses are all less than 1.5 Mq. 

On the other hand, within the nuclear many-body framework, Jiang et al. have shown that 
the inclusion of ring diagrams in the DBHF calculations softens the equation of state of symmetric 
nuclear matter. It would be of interest to extend this to asymmetric nuclear matter, and to see 
to what entent this improvement in the nuclear many-body approach can improve the agreement 
with the astrophysics observations. 
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5 Summary 



In summary, we studied the properties of asymmetric nuclear matter in the formalism of the Dirac- 
Brueckner approach with the Bonn one-boson-exchange nucleon-nucleon interaction. The symmetry 
energy coefficient at the saturation density obtained in this work is about 30 MeV. This is in good 
agreement with the empirical value of about 34±4 MeV, and in agreement with other approaches 



such as the BHF [26] and the variational [25] calculations. The symmetry energy in our study was 
found to increase almost linearly with the density and agrees with the linear parameterization of 
Prakash, Ainsworth and Lattimer ||. At higher densities, the symmetry energy in our calculation 
is considerably larger than those in the BHF and variational calculations. The difference can be 
understood as coming from the both the relativistic effects in the 'kinetic energy' contribution, and 
a strong p-meson coupling in the Bonn potential that increases the 'potential energy' contribution, 
to the symmetry energy. 

We have also applied the resulting equation of state of neutron-rich matter to calculate the 
maximum mass of nucleon star, and we find it to be about 2.1Mq. The corresponding radius 
and central density are R = 10.8 km and p c = 1.08 fm~ 3 , respectively. This maximum mass 
is substantially greater than that measured in any neutron stars. It would be brought down 



considerably if kaon condensation were included in the EOS [39, 41]. We believe our maximum 
mass of 2.1 Mq, without inclusion of kaon condensation, to be the most direct determination of 
this important parameter, since we begin from a microscopic two-body interaction and use the best 
presently available technology to calculate the mass. 
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Figure 1: Left panel: neutron and proton scalar and vector potentials as a function of asymmetry 
parameter, for three densities. Right panel: equation of state of nuclear matter for a number of 
asymmetry parameters. 
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Figure 2: Left panel: energy difference AE as a function of a 2 for several densities with u = p/p 
and po = 0.166 fm~ 3 . Right panel: density dependence of symmetry parameters. 



13 




Figure 3: Left panel: comparison of our results with those of Refs. |25| and ^6[. Right panel: 
comparisons of our results with phenomenological parameterizations of Prakash, Ainsworth and 
Lattimer |H|. 
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Figure 4: Left panel: proton fraction and AE. Right panel: nucleon star properties. 
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